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摘要 : 高 海拔 宇宙 线 观 测 站 (Large High Altitude Air Shower Observatory, LHAASO) 主要 
用 于 探索 高 能 宇宙 线 起 源 、 高 能 天 体 演化 以 及 新 物理 前 沿 等 。 其 中 的 大 面积 缘 子 探测 器 (Muon 
Detector, MD) 阵列 ， 具 有 对 原初 伽 蕊 射线 的 高 鉴别 能 力 ， 能 够 完成 超 高 能 佑 蕊 射线 的 高 灵 
敏 度 探测 。 绿 子 探测 器 通过 阳极 和 第 七 打 拿 极 两 路 信号 共同 读 出 的 方式 ， 实 现 1~10 000 个 
比 子 的 大 动态 范围 探测 。 为 了 减 小 在 阳极 信号 饱和 后 的 区 域 ， 直 接 转 换 为 打 拿 极 信 号 进行 测 
量 的 误差 ， 提出 了 将 阳极 (阳极 未 饱和 部 分 ) 和 打 拿 极 ( 阳极 饱和 部 分 ) 两 部 分 信号 进行 拼接 
测量 的 电荷 测量 优化 方案 。 通 过 数值 模拟 得 出 阳极 与 打 拿 极 之 间 的 连接 方式 ， 可 以 采用 实验 
室 标定 的 阳极 、 打 拿 极 关系 进行 连接 。 数 值 模拟 结果 表明 ， 该 方案 在 1~10 000 个 缪 子 的 大 
动态 范围 内 有 效 降 低 了 测量 误差 ， 误 差 小 于 1% ， 缪 子 谱 结构 得 到 明显 改善 。 利 用 全 波形 信 
号 进行 数据 分 析 的 结果 与 模拟 结果 很 好 符合 ， 验 证 了 该 方案 的 有 效 性 。 
关键 词 : 电荷 测量 ; 高 海拔 宇宙 \ 线 观测 站 ; STRMB; 非 线 性 
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高 海拔 宇宙 线 观 测 站 通过 探测 需 探 测 初级 宇宙 线 粒 子 经 过 广 延 大 气 篮 射 (FExtensive Air Shower, 
EAS) 后 产生 的 次 级 粒子 ， 通 过 重建 、 分 析 给 出 原初 宇宙 线 的 相关 信息 。 地 面 复 射 阵列 (KM2A ) 作为 
高 海拔 宇宙 线 观 测 站 的 主体 阵列 ， 由 电磁 粒子 探测 需 (Electron Detector, ED ) KEJ AZ T IRW 4% z7] 
组 成 ， 用 来 探测 广 延 大 气 簇 射 中 的 次 级 电磁 粒子 和 比 子 ， 从 而 开展 超 高 能 宇宙 线 起 源 和 宇宙 线 物理 的 
研究 。 其 中 ， 缪 子 探测 器 的 布置 总 面积 达 40 000 m°, JEFF 1 188 个 探测 单元 ， 采 用 反射 式 水 切 伦 科 夫 
探测 技术 对 广 延 大 气 徐 射 中 的 纽 子 成 分 进行 探测 。 纪 子 探测 絮 的 大 面积 铺设 极 大 提高 了 对 原初 伽 马 射 
线 的 鉴别 能 力 ， 实 现 了 对 超 高 能 伽 马 射线 的 高 灵敏 度 探测 ， 直 接 有 瞄准 解决 超 高 能 宇宙 线 起 源 这 一 目 
标 。 广 延 大 气 篮 射 次 级 纱 子 含量 测量 有 助 于 我 们 进行 原初 军 宙 线 成 分 区 分 和 高 能 强 相 互 作用 模型 的 鉴 
别 ， 使 得 LHAASO 可 以 精确 测量 “ 膝 ” 区 成 分 能 谱 ， 进 而 破解 “ 膝 ” 区 成 因 这 一 难题 。 缪 子 探测 
需 单 体 设计 如 图 1， 饶 体 上 面 履 盖 2.5 m 厚 的 土屋 ， 吸 收 缪 子 探测 带 中 的 电磁 成 份 ， 提 高 探测 带 探 测 
到 比 子 的 纯度 。 在 圆柱 形 水 泥 色 中 放置 一 个 密封 性 能 良好 且 有 高 反 冉 率 内 层 的 水 从， 水 供 中 灌 入 超 纯 
水 。 纪 子 穿 过 水 袋 时 在 水 体 中 产生 切 伦 科 夫 光 ， 这 些 光子 在 水 中 传播 时 被 内 袋 反 射 ， 最 终 有 一 部 分 光子 
由 安装 在 袋子 上 表面 中 央 的 光电 倍增 管 (Photo Multiplier Tube, PMT) 接收 。 光 电 售 增 管 输出 阳极 和 
打 拿 极 信 号 ， 通 过 同 轴 电 绕 接 入 电子 学 ， 经 过 前 端 模拟 电路 整形 后 送 入 flash ADC 进行 采样 ， 采样 数据 
送 入 现场 可 编程 门 阵列 ( Field Programmable Gate Array, FPGA ) 中 进行 实时 人 处理， 最 终 通过 光纤 输出 le 

23 RM at BOR IA EE 1 ~ 10000 个 粒子 ,为 实现 大 动态 范围 测量 ,采用 阳极 和 第 七 
打 拿 极 共同 读 出 的 方式 ， 阳 极为 1~166 个 粒子 ， 打 拿 极 为 58~10 000 个 粒子 “。 由 于 采用 阳极 加 打 
拿 极 共同 读 出 的 方式 ， 需 要 建立 较为 准确 的 阳极 和 打 拿 极 连接 关系 。 同 时 ， 阳 极 信 号 饱和 后 将 强制 转 
换 为 打 拿 极 信号 ， 此 时 打 拿 极 信号 较 小 ， 存 在 电子 学 非 线 性 较 大 、 测 量 误差 较 大 的 问题 。 为 了 减 小 测 
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量 误差 ， 本 文 提 出 了 一 个 比 子 探测 副 电 答 测 量 优 化 方案 ， 并 利用 数值 模拟 和 全 波形 数据 分 析 其 可 行 
性 ， 验 证 其 有 效 性 。 


Ø13.9 m 


Electronic 
system 


Soil 


PMT 


Anode | Frontent 
analog ===» FADC 
circuit 


J[qvd [BIXBOD 


Frontent 


analog ===» FADC ==} 
Dynode| circuit 


Ultra-pure water 


Concrete cistern 


Electronic system 


Al Be rR MI B24 RSA (LHAASO MD) 
Fig. 1 Schematic diagram of MD unit detector (LHAASO MD) 


1 ZTM Em ee TT 


1.1 电子 学 输入 信号 特征 
目前 观测 到 的 宇宙 线 能 量 跨 10 多 个 数量 级 ， 流 量 强度 跨 30 个 数量 级 ， 它 的 能 谱 总 体 上 服从 简单 
的 备 率 分 布 ， 高 能 宇宙 线 经 过 广 延 大 气 簇 射 形成 的 缪 子 流量 强度 也 基本 呈 备 率 分 布 ”。 
入 则 缘 子 探测 右 的 信号 经 过 光电 倍增 管 接收 转换 为 电信 号 ， 青 输入 电子 学 系统 进行 处 理 。 探 测 带 
的 能 量 分 辨识 反映 了 线性 探测 带 区 分 被 探测 粒子 
能 量 的 能 力 。 图 2 为 单个 比 子 的 积分 电 奏 分 布 情 
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表示 ， 其 中 ，A 为 水 吸收 长 度 ; fA SEE, /为 光子 在 两 次 与 内 袋 碰 撞 的 平均 径 迹 长 度 ; AIG 
在 真空 中 的 速度 ; 为 光 在 水 中 的 反射 率 。 在 实际 情况 中 ， 这 些 参数 基本 稳定 ， 因 此 衰减 时 间 也 基本 
稳定 。 波 形 前 治 上 升 时 间 较 短 (4~9 ns) ， 且 随 缪 子 数 变化 不 大 “。 图 3 为 采用 4 000 SEB BOE AF LY 
化 的 波形 图 ， 前 沿 上 升 较 快 ， 而 后 沿 为 e 指数 衰减 曲线 。 

对 于 平均 波形 ， 其 后 沿 为 较 好 的 。 指数 娶 减 曲线 ， 但 对 于 单个 波形 来 说 ， 光 子 经 过 反射 后 一 波 一 
波 到 达 光 电 售 增 管 并 被 吸收 ， 因 此 ， 后 沿 表现 为 在 平均 波形 附近 晃动 ， 即 波形 特征 晃动 ， 晃 动 幅 度 与 
平均 波形 的 幅度 相关 。 对 2 000 个 波形 的 平均 波形 进行 拟 合 得 到 平均 值 ， 图 4 为 单个 波形 晃动 幅度 与 
平均 幅度 之 比 的 分 布 情况 ， 基 本 服从 高 斯 分 布 。 
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Fig.3 Mean wave form of a single muon Fig.4 Characteristic sloshing distribution of waveform 


1.2 电荷 测量 原理 及 电子 学 特征 

缪 子 探测 需 读 出 电子 学 作为 前 端 电子 学 ， 主 要 用 于 接收 缪 子 探测 需 的 光电 倍增 管 输出 信号 ， 完 成 
电荷 量 O、 时 间 量 了 的 测量 ， 采 样 频率 为 S00 MHz。 同 时 ,每 1s 给 出 一 个 信号 完整 的 波形 供 物 理 分 
析 使 用 。 

电子 学 的 阳极 和 打 拿 极 信 号 采用 两 路 独立 模 数 转换 器 分 别 进行 波形 数字 化 ， 数 字 化 的 信号 采用 数 
值 积分 方法 测量 电荷 量 ， 积 分 窗口 约 600 ns, WEMATZ) 100 ns, LEV 500 ns。 模 数 转 换 器 按时 钟 
节拍 对 成 形 后 的 信号 不 间断 地 进行 瞬间 取样 和 数字 化 ， 得 到 一 系列 数字 化 取样 值 Do, Di, Doi, o, 
D、， 当 触发 判 选 有 效 时 ， 再 对 这 些 数字 化 的 结果 进行 积分 ， 积 分 面积 与 信号 电荷 量 成 正比 ， 得 到 积 
分 面积 即 可 计算 电荷 量 … ” 。 图 5 为 电荷 积分 示意 图 。 
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Fig.5 Charge integration diagram 


由 于 受 模 数 转换 器 的 固有 噪声、 工艺 上 的 缺陷 和 各 种 干扰 噪声 的 影响 ， 电 子 学 本 身 存在 噪声 。 而 
电子 学 台阶 的 性 能 反映 了 电子 学 噪声 水 平 "' ， 表 现 为 电子 学 无 信号 输入 时 ， 输 出 信号 存在 晃动 现象 ， 
如 图 6 为 扣除 电子 学 人 台阶 后 的 晃动 情况 。 

1.3 阳极 与 打 拿 极 之 间 的 转换 关系 (AD ratio) 标定 

对 于 电子 学 最 终 和 输出 的 电荷 量 ， 阳 极 信号 可 以 根据 单 技 电荷 标定 建立 电荷 量 与 纱 子 数 的 转换 关 
系 ， 但 对 于 打 拿 极 信号 ， 则 需要 先 将 打 拿 极 信号 转换 为 等 效 阳 极 信号 后 才能 进一步 转换 为 缘 子 数 ， 因 
K, 我们 需要 建立 阳极 与 打 拿 极 之 间 的 转换 关系 ， 即 AD ratio, AD ratio 分 为 真实 值 与 标定 值 ， 其 中 ， 
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真实 值 为 实验 室 对 光电 倍增 管 进行 标定 得 出 ， 标 定时 采用 的 测量 电子 学 非 线性 及 误差 基本 可 以 忽略 ， 
因此 真实 值 反映 光电 倍增 管 实际 的 AD ratio; 标定 值 则 是 采用 光电 倍增 管 配套 的 谈 出 电子 学 系统 读 出 
数据 进行 标定 ， 具 体 标 定 方法 有 以 下 两 种 : (1) 积分 电荷 (Q) 标 定 法 。 采 用 缪 子 探测 需 输 出 的 积分 电 
谷 进 行 标定 ， 通 过 对 未 饱和 阳极 信号 且 信和 号 较 大 的 阳极 和 打 拿 极 电荷 量 之 比 进行 拟 合 可 以 得 到 电荷 
AD ratio, WHAI7, FP REAR iE A BARA, DAR AIT ERB PM, FH PREY 608 ns, 
在 满足 一 定数 据 量 时 ， 扣 除 电子 学 台阶 后 的 晃动 平均 值 为 0，AD ratio 基本 只 受 电 子 学 非 线 性 影响 ， 
因此 能 真实 反映 阳极 饱和 后 区 域 的 输出 阳极 信号 与 打 拿 极 信号 之 间 的 关系 。(2) 波形 峰值 (peak ) 标定 
法 。 采 用 比 子 探测 癌 输 出 的 全 波形 信号 进行 标定 ,通过 对 未 饱和 阳极 信号 且 信 号 较 大 的 两 路 信号 波形 
进行 寻 峰 ， 得 到 的 峰值 进行 拟 合 可 以 得 到 峰值 AD ratio, 587 AD ratio 的 拟 合 方式 一 致 。 由 于 采用 
单 点 标定 ， 受 电子 学 台阶 园 动 影响 ,波形 的 各 个 点 受 晃动 影响 偏 大 或 仿 小 ， 打 拿 极 寻 到 的 峰 并 不 一 定 
为 真实 峰值 ， 因 此 ， 进 行 标定 后 的 AD ratio 并 不 能 真实 反映 阳极 与 打 拿 极 之 间 的 关系 。 
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Table 1 Basic model parameters 

2.1 模型 构建 及 初 步 验证 Parameter Value 

模型 的 构建 主要 包括 对 电子 学 输入 信号 和 电 Exponents of the muon power law 2. 47 
子 学 的 真实 模拟 ， 通 过 上 述 分 析 确 定 模拟 中 采用 spectrum 
的 基本 参数 如 表 1。 Energy resolution of a single Muse 25% 

采用 转换 方案 测量 时 , PA 极 饱和 后 , 在 信号 The attenuation index of the back edge 200 
较 小 时 AD ratio 采用 标定 值 误 差 较 小 ， 信 号 较 大 of the waveform 
时 采用 真实 值 误差 较 小 两 个 AD ratio 的 切换 点 The growth index at the front of the 4 
与 标定 值 之 间 存 在 简单 的 线性 关系 ， 通 过 标定 值 人 
可 以 确 ce AD ratio 的 HJ 换 点 5 Waveform characteristic sloshing Guas 
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结构 出 现 异 常 如 图 8 (a) hy BEB HF PEG Charge integration window, ns 608 
况 ， 阳 极 饱 和 后 到 AD ratio 切换 点 的 区 间 段 ( 约 DNL Guas (mean, sigma), LSB Guas 

e m [11] 

170~270 FF) 由 于 AD ratio 采用 标定 值 ， 而 We 
标定 值 实际 偏 大 l He J BATAK , 出 现 了 Electronics step shaking Cuas Guas 


1% [XN ik YP: J 、 mean, sigma), LSB ,1.2 
图 中 该 区 间 段 谱 线 下 移 的 现象 ; 300 TA T ME mean Sums) aT (0, 1. 20) 


a 


262 H 


K X 研 与 Ñ 20 48 
为 AD ratio 切换 点 ， 切 换 后 AD ratio RAB, MAXXE, Te AS ev), h Ale 
300 TAT Bi HAZE Ie] EA, BUTS EIA FA 8(b)， REUSE TNA t 


数 十 分 吻合 (<0. 16%) ， 初 步 验证 了 模型 的 可 靠 性 。 
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Fig. 8 The muon power law spectral structure of the conversion scheme 


2.2 AD ratio 标定 分 析 

对 于 波形 峰值 和 积分 电荷 两 种 AD ratio 标定 法 ， 我 们 通过 模拟 和 数据 分 析 了 标定 值 与 标定 区 域 比 
子 数 之 间 的 变化 关系 。 

(1) 图 9(a) 为 AD ratio 采用 波形 峰值 标定 法 进行 标定 时 ，AD ratio( 按 真实 值 进行 归 一 化 ) bi 
数 的 变化 情况 。 比 子 数 较 小 时 ， 打 拿 极 峰值 受 电子 学 台阶 园 动 的 影响 相对 较 大 ， 寻 峰 后 的 打 拿 极 总 体 
峰值 偏 大 ， 拟 合 的 AD ratio 偏 小 。 随 着 刻 子 数 增 大 ， 受 电子 学 台阶 时 动 的 影响 相对 减 小 ，AD ratio 逐 
渐 增 大 。 尽 管 在 大 信号 时 采用 峰值 标定 的 AD ratio 接近 真实 值 ， 但 仅仅 是 电子 学 台阶 晃动 造成 的 巧 
合 ， 当 晃动 幅度 改变 时 ， 标 定 值 也 发 生 改 变 ， 因 此 并 不 能 真实 反映 阳极 和 打 拿 极 的 关系 。 
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图 9 AD ratio 标定 模拟 及 数据 


Fig.9 AD ratio calibration simulation and data 
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(2) Al 9(b) 4 AD ratio 采用 积分 电 丛 标定 法 进行 标定 时 ，AD ratio EZTA Eki, BR 
较 小 时 ， 打 拿 极 受 非 线性 影响 较 大 ， 积 分 电荷 总 体 俩 小 ， 因 此 拟 合 后 AD ratio WAK, tia BF By 
大 ， 非 线性 逐渐 减 小 ，AD ratio 也 逐渐 减 小 。 由 于 电子 学 最 终 输 出 值 为 电 谷 量 ， 此 时 利用 电荷 量 进行 
标定 的 AD ratio 正好 可 以 真实 反映 该 人 射 舱 子 数 附近 阳极 与 打 拿 极 之 间 的 关系 。 模 拟 与 数据 之 间 存 在 
一 定 偏 差 .这 是 由 于 打 拿 极 信号 波形 中 包含 大 量 小 幅度 信号 ， 由 电子 学 微分 非 线性 (Differential 
nonlinearity, DNL) 可 以 得 出 ， 此 时 非 线 性 变化 很 大 ， 因 此 相应 的 AD ratio 随 纱 子 数 的 变化 趋势 也 存在 
差异 。 如 采 要 准确 测量 电荷 ， 应 该 尽 可 能 避免 采用 这 部 分 小 幅度 的 信号 ， 这 也 是 采用 拼接 方案 进行 电 
集 测 量 的 一 个 原因 。 
2.3 ” 电 三 测量 拼接 方案 及 模拟 

现 有 的 电荷 测量 方案 ， 即 转换 方案 ， 电 子 学 最 终 输 出 为 阳极 和 打 拿 极 积分 电 和 人 向， 在 阳极 信号 饱和 
后 ， 只 能 采用 打 拿 极 积分 电 丛 进 行 计算 ， 在 阳极 饱和 区 域 ， 打 拿 极 信号 较 小 ， 受 电子 学 非 线 性 影响 ， 
打 拿 极 积分 电 三 误差 较 大 。 因 此 ， 通 过 分 析 提 出 采用 阳极 未 饱和 部 分 的 信号 和 阳极 饱和 部 分 对 应 的 打 
拿 极 信号 ， 两 部 分 信号 拼接 进行 测量 ， 即 拼接 方案 。 如 图 10, 饱和 点 以 下 的 部 分 采用 阳极 信号 计算 
电 俩 量 ， 而 饱和 点 以 上 的 部 分 只 有 打 拿 极 可 以 测 得 ， 这 部 分 采用 打 拿 极 信号 计算 电 全 量 。 由 于 采用 两 
部 分 信号 测量 ， 阳 极 部 分 的 信号 非 线 性 有 影响 很 小 ， 打 拿 极 部 分 的 信号 采用 波形 信号 较 大 的 部 分 ， 非 线 
性 影响 也 较 小 ， 可 以 有 效 降低 总 体 测 量 误差 。 
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图 10 拼接 方案 示意 图 
Fig. 10 Splicing scheme Diagram 


图 11 RADE DT RET OBE, y 
阳极 饱和 区 域 的 谱 结 构 得 到 明显 改善 。 由 于 拼接 
方案 采用 未 饱和 阳极 信号 和 多 和 阳极 信号 对 应 的 


打 拿 极 信号 ， 该 部 分 信号 非 线 性 较 小 (<3%)， p” sonst 
因此 可 以 直接 采用 真实 值 进行 计算 。 图 12 为 两 。 
种 方案 的 电荷 积分 测量 误差 模拟 ， 转 换 方案 的 10 

AD ratio 在 缪 子 数 较 小 时 采用 标定 值 ， 由 2.2 节 $ 

图 9(b) 可 以 得 出 ， 随 着 织 子 数 增加 ， 应 该 采用 10 


的 AD ratio 逐渐 降低 至 真实 值 ， 因 此 使 用 标定 值 

计算 的 积分 电 集 误差 逐渐 增 大 ， AD ratio FZT | 

数 较 大 时 采用 真实 值 此 时 计算 后 的 积分 电荷 实 70 8090100 200 300 400 500 600 1 000 
Number of 

际 偏 小 ， 但 随 着 缪 子 数 的 增 大 误差 逐渐 降低 。 拼 Die dee ee 

接 方案 则 只 采用 真实 值 计 算 ， 此 时 积分 电荷 总 体 nop BARA hn E o S 

偏 小 ， 随 着 缪 子 数 增 大 ， 打 拿 极 部 分 的 比例 不 断 a d P 

上 升 ， 误 差 逐渐 增 大 ， 但 非 线性 的 不 断 降 低 使 得 

误差 到 达 最 大 后 又 逐渐 降低 。 总 之 ， 拼 接 方案 在 阳极 饱和 区 域 的 误差 得 到 有 效 控制 ， 误 差 小 于 1%. 
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3 BP Pel te Ba Sore 


ZTR ai BE a US SC EO, 
括 阳 极 信号 和 打 拿 极 信号 ， 波 形 窗口 为 5 ps, fh 
发 前 约 500 ns， 触 发 后 约 4500 ns， 利 用 这 些 信 
号 可 以 完成 不 同 目标 的 物理 分 析 ， 图 13(a) 为 单 
次 事例 的 阳极 波形 图 ， 图 形 上 方 平整 部 分 为 阳极 
饱和 区 域 ， 对 应 的 打 拿 极 波形 如 图 13(b)。 采 用 
拼接 方案 进行 电荷 量 计 算 ， 电 荷 积 分 窗口 608 ns, 
过 国 前 100 ns, WW a 508 ns, AD ratio 使 用 真 
KE, FENN Fie uR 14， 谱 结构 得 到 明显 
改善 ， 同 时 ， 谱 指数 与 数值 模拟 情况 十 分 吻合 
(<0.24% ) ， 验 证 了 该 方案 的 可 行 性 。 
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Fig. 13 Single case waveforms at anode saturation 


电荷 测量 方案 的 模拟 及 数据 分 析 如 表 2， 两 
种 方案 数据 与 模拟 的 谱 指数 十 分 吻合 ， 拼 接 方案 
与 原初 模拟 谱 更 接近 ( <0. 29% ) 。 数 据 结 果 中 的 
X 比 模拟 情况 更 大 ， 这 是 由 于 使 用 了 1 000 个 探 
测 需 近 1 年 的 数据 进行 分 析 ， 不 同 探测 器 间 存 在 
一 定 差异 ， 同 一 探测 器 不 同 环境 条 件 下 也 存在 一 
定 差异 ， 但 在 误差 范围 内 ， 模 拟 与 数据 结果 十 分 
吻合 。 采 用 拼接 方案 时 可 直接 使 用 真实 值 作 为 
AD ratio， 其 电荷 测量 误差 <1% ， 此 时 标定 值 只 
作为 判断 探测 器 是 否 存在 异常 的 一 个 指标 。 
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3 期 
表 2 电荷 测量 方案 对 比 
Table 2 Comparison of charge measurement schemes 
Measurement scheme Spectral index x AD ratio Error/% 

Conversion scheme (simulation ) 2. 484 1. 681 Pee ied thon value zÁ 

Conversion scheme ( data) 2. 488 2. 426 
Splicing scheme (simulation ) 2.477 1. 129 True value “4 

Splicing scheme ( data) 2. 471 1. 408 


局 灵敏 度 探测 的 基础 。 针 对 阳极 、 打 拿 极 连接 问题 ， 本 文通 过 理论 分 析 和 数值 模拟 相 结 合 的 方式 ， 得 
出 拼接 方案 可 以 直接 采用 实验 室 标定 的 AD ratio 真实 值 进行 连接 ， 探 测 带 反馈 数据 标定 的 AD ratio 标 
定 值 只 作为 判断 探测 带 是 否 存在 异常 的 一 个 指标 ， 有 效 避 免 了 AD ratio 标定 及 切换 可 能 引入 的 误差 ， 
简化 了 电荷 积分 的 计算 过 程 。 通 过 数值 模拟 ， 我 们 分 析 验 证 了 拼接 方案 可 以 有 效 降 低 阳 极 饱和 区 域 电 


fy Wei, SAUER ARTA AIL, IRA RK 4% 降 低 至 最 大 1%， 实 现 了 在 1~ 10 000 BF 
的 大 动态 范围 的 精确 测量 。 同 时 对 纪 子 探测 器 反馈 的 全 波形 信号 采用 拼接 方案 计算 电荷 量 ， 织 子 谱 结 
构 和 谱 指 数 与 模拟 情况 十 分 吻合 ， 验 证 了 该 方案 的 有 效 性 。 
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Optimization of Charge Measurement for LHAASO Muon Detector 
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2. Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 100049, China, Email: zuoxiong@ihep.ac.cn) 


Abstract: The Large High Altitude Air Shower Observatory (LHAASO) is mainly used for exploring the 
origin of high-energy cosmic rays, the evolution of high-energy celestial bodies and the research of new 
frontiers of physics. The large area Muon Detector (MD) array has high identification ability of primary gamma 
rays and can detect ultra-high energy gamma rays with high sensitivity. MD achieves a large dynamic range 
detection of 1-10 000 muons through the common readout of the anode and the seventh dynode. In order to 
reduce the error of directly converting the anode signal into the dynode signal in the saturated area, an 
optimization scheme of charge measurement was proposed by combining the anode (the unsaturated part of the 
anode ) and the dynode (the saturated part of the anode) signals. The connection mode between the anode and 
the dynode is obtained by numerical simulation, and the relationship between the anode and dynode calibrated 
in the laboratory could be used for connection. Numerical simulation results show that the proposed scheme can 
effectively reduce the measurement error in the large dynamic range of 1-10 000 muons, the error is less than 
1%, and the muons spectrum structure is obviously improved. The results of data analysis with full waveform 
signal are in good agreement with the simulation results, which verifies the effectiveness of the scheme. 


Key words: charge measurement; LHAASO; MD; nonlinear 


